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We present measurements of silica optomechanical resonators, known as bottle resonators, passively cooled 
in a cryogenic environment. These devices possess a suite of properties that make them advantageous for prepa- 
ration and measurement in the mechanical ground state, including high mechanical frequency, high optical and 
mechanical quality factors, and optomechanical sideband resolution. Performing thermometry of the mechanical 
motion, we find that the optical and mechanical modes demonstrate quantitatively similar laser-induced heating, 
limiting the lowest average phonon occupation observed to just ~1500. Thermalization to the 9 mK thermal 
bath would facilitate quantum measurements on these promising nanogram-scale mechanical resonators. 


I. INTRODUCTION 


Cavity optomechanical systems, consisting of high-quality 
optical cavities coupled to mechanical resonators, offer a 
promising route to making precision measurements in both the 
applied and fundamental science domains [1]. In dispersively- 
coupled systems, the motion of the mechanical resonator 
shifts the resonance frequency of the optical cavity, which 
can in turn be observed as periodic changes in the ampli- 
tude or phase of light traversing the cavity. Optomechanics 
thus provides an extremely sensitive readout for micro- and 
nano-mechanical resonators, enabling their use as exquisite 
sensors of a variety of phenomena on small scales, including 
displacement [2], force [35], torque [6] [7] and acceleration 
[8]. This readout sensitivity has also motivated fundamen- 
tal searches for quantum properties of nanomechanical res- 
onators [9] at, or near, their vibrational ground state where 
mechanical motion originates from quantum zero-point fluc- 
tuations. The hybrid nature of optomechanical systems ad- 
ditionally makes them desirable for applications in quantum 
information processing architectures [10], and as a resource 
for entangled [11], squeezed or other non-classical states 
[13). 


Noise from thermally-driven oscillations of mechanical res- 
onators poses a significant obstacle to performing many of 
these proposed experiments. Mechanical resonators with fre- 
quencies in the kHz to GHz range have thermal occupancies of 
billions to thousands of phonons at room temperature, which 
can easily drown out any quantum signature. Reducing this 
thermal occupation is thus of great importance. The optome- 
chanical interaction can be exploited to actively cool the me- 
chanical mode into or near its ground state but a fun- 
damental limit to this process is imposed by the temperature of 
the mechanical resonator’s bath, necessitating cryogenic pre- 
cooling [I]. Active optomechanical cooling furthermore re- 
duces the quality factor (Qm) of the mechanical resonator, de- 
creasing the signal-to-noise ratio at low temperatures [I]. For 
these reasons, we focus on passively cooling our resonators, 
which additionally facilitates the use of sensitive optomechan- 
ical systems for probing low temperature phenomena such as 
superfluidity and superconductivity [19]. 

In previous experiments, optomechanical systems have 


been passively cooled using both helium flow cryostats 
and dilution refrigerators [11] [21]. Our system uses a di- 


lution refrigerator for its lower achievable base temperature 
and, in contrast to the system presented in Ref. uses the 
highly-efficient tapered-fiber coupling method [22]/23], which 
is compatible with both bulk and on-chip optomechanical sys- 


tems [24] [25]. 


Here, we present measurements of so-called silica bottle 
resonators [26] /27], shown in Figure[I] which have ellipsoidal 
shapes and exhibit optical whispering gallery modes (WGM) 
in the infrared and visible range, a simulation of which is 
shown in Figure [I{b). They are interesting structures for their 
high optical quality factors (Qo up to 108) and the tunability 
of their mode structure [27]. Additionally, these bottle res- 
onators have mechanical breathing modes with frequencies in 
the 50 to 250 MHz range, Figure [I{c), corresponding to av- 
erage phonon occupancies of just 3.3 and 0.4, respectively, 
at a temperature of 9 mK (the base temperature of our dilu- 
tion refrigerator). The significant modal overlap of the optical 
and mechanical modes, which are both localized at the equa- 
tor, leads to a large optomechanical coupling, on the order of 
G/2n ~ 10 GHz/nm. Furthermore, the combination of high 
optical quality factor and high mechanical frequency places 
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FIG. 1: (a) Photograph of a string of bottle resonators along with 
the tapered fiber for optomechanical measurement. Scale bar is 100 
jum. (b) Simulation of a whispering gallery mode optical resonance, 
showing a cross-section of the bottle with the optical field localized 
along the equator. (c) Representative mechanical spectrum, mea- 
sured at room temperature. The three most prominent resonances 
are indicated by the corresponding mechanical mode simulation. 


these devices in the sideband-resolved regime, enabling many 
of the important tools of optomechanics [I]. 


Il. EXPERIMENTAL DETAILS 


We fabricate our bottle resonators by applying tension to 
each end of a single mode optical fiber while melting it with 
a COz laser. This process results in a string of resonators, as 
shown in Figure[I{a), separated from each other by thin stems. 
By controlling the intensity of the CO2 laser and the length of 
the pull time, we create resonators with symmetric shapes and 
extremely thin supporting stems. This reduces phonon tun- 
nelling through the stems and increases the mechanical qual- 
ity factor [28]. 

Tapered fibers are used to couple light into the 
WGMs of the bottles. These fibers have a small core (on the 
order of the wavelength of the light) and no cladding, save 
the medium surrounding the fiber. Our tapers are created by 
heating an optical fiber with a hydrogen torch and pulling it 
until a minimum diameter of ~1 um is reached [25], result- 
ing in a large evanescent field. This allows efficient coupling 
of light to WGM resonators through frustrated total internal 
reflection. 

We passively cool the taper and bottle resonators using a 
dilution refrigerator with a base temperature of 9 mK [24]. To 
facilitate taper-resonator coupling at cryogenic temperatures, 
we have built an optomechanical coupling system within the 
inner vacuum can (IVC), on the mixing chamber plate, as pic- 
tured in Figure[2[a). The bottles are mounted on top of a stack 
of nanopositioning stages, which allow full three-dimensional 
control over the position of the resonators. They are thermally 
anchored to the mixing chamber through a series of oxygen- 
free high-conductivity copper plates and braids. The taper 
is epoxied and mounted to a large Invar block, which min- 
imizes the effects of thermal contractions. A homebuilt low- 
temperature microscope enables real-time in situ imaging 
of the taper-resonator system with a resolution of ~1 jum (Fig- 
ure [2{b)). We inject light from a tunable diode laser housed at 
room temperature into the tapered fiber and detect the trans- 
mitted light with a fast photodetector. The optical and me- 
chanical properties of the bottles can thus be measured from 
the low- and high-frequency parts of the photodetector volt- 
age, respectively. 


Il. OPTICAL PROPERTIES 


Characterization of the near-infrared optical modes is per- 
formed by scanning a tunable diode laser (1500-1630 nm) 
with angular frequency @,, across a resonance, resulting in the 
following expression for the low-frequency part of the trans- 
Mission Aoyt (in units of photons per second) through the ta- 
pered fiber 


Aout = Nin — KONcav- (1) 


FIG. 2: (a) Optomechanical coupling apparatus on the mixing cham- 
ber stage of our dilution refrigerator. (b) Image of a bottle resonator 
and tapered fiber taken with the low-temperature microscope. Scale 
bar is 100 um. 


Here fin = Pi, /h@, is the rate of photons passing through 
the tapered fiber, for injected laser power P,, and reduced 
Planck constant i. The linewidth of the resonance is given 
by K = Ko + Kex where Ko is the intrinsic decay rate and Kx 
is the rate at which photons are exchanged between the bot- 
tle and tapered fiber. Our setup allows us to sensitively con- 
trol the relative strengths of these intrinsic and extrinsic decay 
rates [24]. For the measurements that follow, the resonator is 
operated in the slightly-undercoupled regime (Kex S Ko). 
The number of photons inside the cavity, 


Nin Kex 
Ncay = Poe y, (2) 
+A 
is a function of the laser detuning A = @, — @, from the cavity 
resonance, @p). For the whispering gallery modes in a bottle 
of radius R, the angular optical resonance frequencies @, are 
given by 


lc 


(3) 


a Ren 
The integer / is a mode label for a particular resonance, c is 
the speed of light in vacuum, n is the refractive index of the 
bottle and c; is a geometric factor which accounts for the fact 
that the optical mode is not perfectly confined to the surface 
of the resonator [29]. 

For relatively low injected optical powers (Pin < 250 nW), 
the optical resonances can be measured without heating the 
mixing chamber thermometer and have Qy of 10°-107, which 
are consistent with their room-temperature quality factors. At 
larger injected powers, light lost from the taper (~30%) heats 
the cryogenic environment [24]. Furthermore, the high Q,, 
and correspondingly long photon lifetime, in the cavity leads 
to an increased absorption of laser photons in the silica. 

This absorption generates a local heating of the silica 
within the optical mode volume, which in turn causes both 
a change in the bottle dimensions through thermal expan- 
sion and a change in the refractive index through the thermo- 
refractive effect. We can then rewrite Equation with 
Taylor-expanded expressions for the radius and the refractive 


index [29], 
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FIG. 3: (a) Simulated normalized taper transmission (black) and 
number of intracavity photons (blue) for a linear optical resonance 
with parameters comparable with those seen in experiment («/27 = 
50 MHz, kex/2m = 20 MHz and Py, = 4 LW). (b) Simulation 
of a nonlinear resonance, with the same parameters as above and 
Aw /2m = —5 x 10~° Hz. All real solutions to Equation (5) are plot- 
ted in black, while the transmission for a forward (reverse) scan of 
the laser is plotted in green (red). 
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where T is the equilibrium temperature of the glass in the ab- 
sence of laser heating, AT is the temperature increase caused 
by optical absorption, and € is the linear thermal expansion 
coefficient. Although higher-order terms become important 
under some conditions [30], here we keep only terms to first 
order in AT and for simplicity write b(T) = ¢++ on 
Incorporating the modified @p into Equation (2), we have 
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(5) 


Ncay = 


where the thermo-optical effects have been written into a non- 
linear detuning parameter A,, (31) [32]. This parameter can be 
thought of as the shift in the optical resonance frequency per 
intracavity photon, and is related to the parameters in Equa- 
tion (4) by Anitcay = b(T)AT, under the assumption that the 
temperature gradient AT is proportional to the number of in- 
tracavity photons. As we will see, the exact constant of pro- 
portionality will depend on the resonator’s heat dissipation. 
Equation (5) is now a nonlinear function of nay, giving rise 
to up to three distinct real solutions for appropriate values of 
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FIG. 4: Transmission measured across a bottle optical resonance 
with a wavelength of 1516 nm and a temperature of 4.2 K. For similar 
injected powers, the resonance exhibits different behavior when the 
bottle is (a) immersed in helium exchange gas and (b) surrounded by 
vacuum. Fits to the data (black line) yield nonlinear detuning param- 
eters A, /2% of approximately (a) —3 x 10-7 Hz and (b) —2 x 10-© 
Hz. 


An. Experimentally, this is manifest as a bistability and hys- 
teresis in the resonance shape, which is dependent on the scan- 
ning direction of the laser (33) [34], as illustrated in Figure 3] 
As the laser is tuned closer to the optical resonance, more pho- 
tons are coupled into the cavity, causing a greater shift in the 
resonance frequency, such that it becomes distorted from the 
linear shape shown in Figure Bla). 


This nonlinear optical behaviour is observed in our exper- 
iments, as shown in Figure [4] for a R= 25 um bottle. Using 
the cubic equation to exactly solve Equation (5), we fit our 
experimental data to extract Aj; under various experimental 
conditions. In particular, we find that injecting a large quan- 
tity of helium exchange gas into the IVC at 4.2 K reduces Ajj 
by an order of magnitude (A, /27 ~ —3 x 10~’ Hz compared 
to —2 x 10-© Hz in vacuum). We furthermore find that A, 
at the fridge base temperature (Ay /27  —3 x 10~° Hz) is 
comparable to that at 4.2 K in vacuum. 


Although it will not be discussed here, it is worth noting 
that at intermediate pressures of helium exchange gas, we ob- 
serve a higher-order thermo-optical nonlinearity. This results 
in up to five real solutions to a modified form of Equation (5), 
and presents itself in the experiment as a multistability in the 
resonance profile. This effect is thought to originate from a re- 
versal of the thermo-refractive effect at low temperatures [30], 
although the nonlinear temperature dependence of the thermal 
expansion coefficient may also contribute. 


IV. MECHANICAL PROPERTIES 


The high-frequency part of the taper transmission encodes 
the mechanical motion of the bottle resonator. We collect this 
signal by recording a time trace of the high-pass filtered pho- 
todetector voltage with a fast analog-to-digital converter. The 
R= 25 um bottle studied has several mechanical resonances, 
with the most prominent being at 55, 85 and 109 MHz, with 
room temperature quality factors of ~ 10+. At low temper- 
atures, we focus on the lowest-frequency (55 MHz) mode, 
which has an effective mass meg = 64 ng and optomechani- 
cal coupling strength G/2a ~ 8 GHz/nm. 

It is important to note that, despite significant efforts to 
thermally anchor the bottle resonators to the base plate of the 
fridge, the thermometers used to measure the temperature of 
the mixing chamber will not provide an accurate measure of 
the mechanical mode temperature. The thermally-insulating 
nature of silica along with the very thin connections between 
the bottles prevent efficient conduction of heat, while incident 
laser light used to measure the mechanics leads to local heat- 
ing of the resonator. Dynamical back-action effects in the op- 
tomechanical interaction can also lead to mode-specific heat- 
ing or cooling |L] which is not reflected in the temperature of 
the bulk silica. 

To independently determine the temperature, we exploit the 
thermally-driven motion of the resonator. In this case, the res- 
onator’s spectral response is quantified by its one-sided power 
spectral density (PSD). For a resonator with effective mass, 
Merf, at temperature T,, the displacement PSD is given by 


Sec() = Jim *|X(o))? 
7 4kpTQm (6) 
mest Om [(Q%, — 07)? + (@Qm/Om)?]’ 


for a sufficiently long measurement time Tt (t >> 27/Q,). 
Here, kg is the Boltzmann constant, Q,, is the angular fre- 
quency of the mechanical resonance and X(@) is the Fourier 
transform of the resonator’s position x(t). It is evident that the 
amplitude of the displacement PSD (4kgT /megr) scales lin- 
early with T and thus gives a direct measurement of the mode 
temperature. 

In our experiments, we measure this displacement response 
by Fourier transforming the photodetector voltage, squaring 
and dividing the result by the measurement bandwidth. The 
voltage PSD is given by the sum of a detection system- 
dependent noise floor S}¥, and the transduction of S,.(@) 
through the optical cavity and photodetector, 


2 
Syv(@) = Spy + fig, J (Qn, A)So( 0). (7) 

m 
The transduction of S,.(@) scales with the square of the rate 
of injected photons, as well as the square of the ratio of 
the optomechanical coupling strength to the mechanical reso- 
nance frequency. All of the detuning dependence of Syy(@) 
is contained in @(Qm,A) (see Appendix), a function which 
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FIG. 5: Schematic of the mechanical mode thermometry technique. 
Modulation of laser light by an electro-optic modulator (EOM) gen- 
erates optical sidebands on the laser spaced by Qmoq. The motion of 
the optomechanical resonator (OMR) imposes additional sidebands, 
at Qm. Both sets of sidebands are detected by the photodiode (PD), 
resulting in a high-pass filtered (HPF) spectrum that has two distinct 
peaks corresponding to the phase modulation (54.5 MHz) and the 
mechanical motion (55 MHz). 


describes the transduction of the resonator’s fluctuations in 
position into fluctuations in first taper transmission and then 
voltage. 

Precise knowledge of the laser detuning, along with the op- 
tomechanical coupling strength G, injected optical power and 
the gains of the detection electronics, would allow the extrac- 
tion of S,.(@), and hence the temperature, from Equation (7). 
This is not always experimentally feasible, so we instead in- 
directly calibrate Syy(@) by injecting an all-optical signal of 
known energy into the system [37]. The general procedure for 
doing so is outlined in Figure [5] An electro-optic modulator 
(EOM) driven at angular frequency Qy,5q phase modulates the 
input laser light, generating optical sidebands on the laser at 
frequencies @, + kKQmoa for integer k. If the phase modula- 
tion depth, B, is small (8 = Vo /Vz < 1 for applied voltage 
Vo and EOM half-wave voltage Vz), only the first order side- 
bands (k = 1) need be considered. If Quyoq is chosen such 
that |Qmod — Qm| >> T where T = Qn/Qm is the resonance 
linewidth, the PSD of the photodetector signal becomes (see 
Appendix) 
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where S4g(@) is the spectrum of the applied phase modula- 
tion. This spectrum features two distinct peaks, as shown in 


Figure 5]at Qmog and Qm. 

If we additionally assume that Qmod & Qm, such that the re- 
sponse of the detection electronics is flat, then &(Qmoa,A) © 
0(Qm,A) and the phase modulation signal and mechanical 
motion are transduced similarly through the optical cavity. 
Although the strength of each peak in Syy(@) depends sensi- 
tively on the injected optical power and the laser detuning, this 
dependence is identical between the two peaks, such that it 
simply provides an overall scale for the spectrum. In contrast, 
any change in T will affect only the mechanical resonance at 
Qm through S,,(@). If the phase modulation conditions are 
held constant, we can then find the temperature of the mode 
through the ratio 


Rp =" eT 6) 
Amod 
where Amoq is the integrated area under the phase modula- 
tion peak and A, is the area under the mechanical resonance. 
For the purpose of these calculations we subtract the detection 
noise floor, ae. 

We drive our EOM at Qmoa/27 = 54.5 MHz and measure 
Rr at liquid helium temperature (4.2 K), where a copious 
quantity of exchange gas in the IVC ensures good thermal- 
ization of the bottle resonator with the outer helium bath. The 
result, R4.2, fixes the temperature measurement scale. The me- 
chanical mode temperature is thus given by 


x 4.2 K. (10) 


Measurements of the bottle resonator made under various 
experimental conditions at low temperatures are shown in Fig- 
ure[6] The laser was scanned across the optical resonance and 
the low- and high-frequency parts of the taper transmission 
were recorded simultaneously. From left to right, measure- 
ments were made at 4.2 K in exchange gas and in vacuum, as 
well as in vacuum with the fridge operating at its base temper- 
ature of 9 mK. In all cases, an injected power of P;, = 24 uW 
was used. 

As the laser is tuned to the center of the optical resonance, 
the low-frequency transmission decreases and more photons 
are coupled into the bottle. There is a corresponding decrease 
in the mechanical resonance frequency, which is accompanied 
by an increase in the mode temperature. In exchange gas, the 
temperature increase is small, on the order of a few kelvin, 
and the maximum relative frequency shift amounts to approx- 
imately —0.04%. In contrast, when the IVC is evacuated, 
we see strikingly similar behavior regardless of whether the 
fridge is operated at 4.2 K or base temperature. In both cases, 
the temperature increases to approximately 40 K, while the 
mechanical frequency decreases by 0.3%. We also observe a 
significant decrease in Qn, from 1100 to 600, as the laser is 
tuned to the center of the optical resonance. 

For a more direct comparison, the mechanical resonance 
frequency and inverse Qm are plotted versus measured me- 
chanical mode temperature in Figure [7| for a number of in- 
jected optical powers. Black triangles indicate measurements 
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FIG. 6: (From top to bottom) The low-frequency transmission of 
light through the tapered fiber traces out the optical resonance pro- 
file (top), while information about the mechanical resonance at 55 
MHz is contained in the power spectral density of the high-frequency 
transmission (plotted here on a log scale). From fits to these spectra, 
we extract the resonance frequency (red), quality factor (green) and 
temperature (cyan) of the mechanical mode across the optical reso- 
nance. 


taken at 77 K (liquid nitrogen) and 295 K (room temperature), 
where the bottles were thermalized using nitrogen exchange 
gas. 


V. DISCUSSION 


Our measurements reveal key information about heat dissi- 
pation in the bottle resonators. Despite efforts to thermally an- 
chor the bottles to the base temperature of the fridge, the low- 
est achieved mode temperature was approximately 4 K, cor- 
responding to an average occupation of just ~ 1500 phonons. 
This temperature was reached when helium exchange gas was 
added to the IVC, indicating that the gaseous helium facili- 
tated convection between the resonator and the bath of liq- 
uid helium. In vacuum, regardless of whether the fridge was 
operated at base temperature or liquid helium temperature, 
the temperature of the bottle was raised upon optomechani- 
cal measurements, as shown in Figures [6] and [7] This is ev- 
idence of the intrinsically weak thermal connection between 
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FIG. 7: (a) Frequency and (b) inverse Qm versus calibrated mode 
temperature. Lines are guides to the eye and black triangles are data 
points at 77 K (liquid nitrogen) and 295 K (room temperature) in 
nitrogen exchange gas. 


the mixing chamber and the silica bottle resonator. 

This heating arises as a result of the absorption of laser light 
into the bulk silica. The degree of nonlinearity in the opti- 
cal resonances, quantified through the parameter Ajj, serves 
as a measurement of the temperature change generated by 
the absorbed light. We find that Aj) is nearly an order of 
magnitude larger in vacuum than in exchange gas, since the 
lack of convection greatly reduces the efficiency with which 
the bottle can dissipate heat. These results are entirely con- 
sistent with measurements of the mechanical frequency shift 
with temperature. As the laser is scanned across an optical 
resonance, we find that the mechanical resonance frequency 
exhibits the same dependence on laser detuning as the taper 
transmission. We furthermore observe a relative frequency 
shift in vacuum that is nearly ten times larger than that in ex- 
change gas. Since the optical and mechanical modes occupy 
nearly the same volume within the bottle structure, heating of 
the optical modes by laser absorption is reflected equally in 
the mechanical mode. 

Finally, we observe a significant dependence of the me- 
chanical quality factor on temperature, with a shape that is 
characteristic of phonon-coupling to configurational two-level 
systems in glass [38]. At high temperatures, thermal activa- 
tion over a potential barrier allows transitions between the two 
configurations; as the temperature is lowered, thermal activa- 
tion ceases and quantum tunneling can occur. At still lower 
temperatures, quantum tunneling is also forbidden and the 
mechanical quality factor increases dramatically. From the 


shape of Figure [7{b), we deduce that our system resides in 
the region of thermal activation. Lower temperatures would 
thus allow us to increase the mechanical quality factor far be- 
yond its room temperature value. There is also an increase in 
the mechanical damping rate when the IVC is filled with ex- 
change gas; we attribute this to interactions with a thin film of 
liquid helium on the surface of the bottle. This may be in fact 
be a promising tool to study ultra-thin superfluid films 
via interactions with the optical and mechanical modes of res- 
onator. We note that the general shapes of both Figure [7{a) 
and (b) are in agreement with the observations of silica toroid 
resonators reported in Ref. 


VI. CONCLUSION 


We have demonstrated passive cooling of a 64 ng optome- 
chanical resonator down to just ~1500 phonons. Further cool- 
ing is prevented by the inability to dissipate the heat caused 
by optical absorption, exacerbated by the high Q, of the sil- 
ica resonator. The scale of this optical absorption was found 
to be in excellent agreement with the degree of heating in the 
mechanical mode. This comparison was enabled by optome- 
chanical mode thermometry, detailed in the Appendix, which 
is now an important tool for quantum optomechanics. Future 
experiments will focus on improving the coupling of these 
sideband-resolved optomechanical resonators to the thermal 
bath. One possibility to achieve this is to use a local reservoir 
of helium as a heat link to the 9 mK dilution refrigerator en- 
vironment, which would reduce optically-induced heating of 
the resonator [40]. In particular, it would be intriguing to use 
liquid helium, which has been shown to provide an excellent 
thermalization medium for micro-electromechanical systems 
(MEMS) down to ~60 mK [41]. Furthermore, the MEMS in 
Ref. 4]|regain their vacuum mechanical dissipation levels at 
T S100 mK, due to the temperature-dependent phonon occu- 
pation in the superfluid state. Successful thermalization of the 
presented nanogram-scale microresonators to 9 mK would re- 
sult in average phonon occupancies of 7 S 3, while maintain- 
ing the high mechanical Q and sideband-resolved nature of 
these optomechanical devices, opening up the door for ground 
state cooling and further quantum optomechanical pro- 
tocols [11]. Finally, extension of these cryogenic silica mi- 
croresonators to doped optical glasses would also enable new 
quantum technologies, such as photonic memories for quan- 
tum cryptographic networks [42]. 
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VII. APPENDIX: MECHANICAL MODE 


THERMOMETRY 
A. Solution of the Optomechanical Cavity 


The equation of motion for the intracavity optical field a(t) 
in an optical cavity that is coupled to a mechanical resonator 
with strength G = d@,/dx is given by 


ios (5 a iGx(t)) a(t)+.Kexsin(t), (CD 


in a frame rotating at the laser frequency @,. Here, sin(t) is 


the input optical field, normalized such that |sin(t)|* = Ain = 
P,,/h@. The field output by the optical cavity is then 


Sout(t) = Sin(t) — /Kexa(t). (12) 


This system can be solved in the stationary regime (see for 
example, Ref. |1) but since we intend to inject a time-varying 
input field, we solve Equation without assuming a sta- 
tionary state. As in Ref.{14] we assume that a(t) can be written 
as the sum of the solution ay (t) to the associated homogeneous 
problem and a particular solution a,(t), 


a(t) = ap(t) +ap(t). (13) 


We further make the assumption that the particular solution 
can be written as ap(t) = ay(t) f(t) where f(t) is a yet-to-be- 
determined function of time. 

We begin by solving for a,(t) by taking s;,(t) = 0. The 
resulting differential equation is 


ats (5 ~iA- iGx(t)) ay(t), (14) 


with the solution 


an(t) = apexp |- é = iA) t +16 [x(oat| (15) 


where ao is an amplitude set by the initial conditions of the 
problem. We choose the form 


x(t) = xge!"/? cos Qmt (16) 


for a damped harmonic oscillator with angular frequency Qy, 
damping rate I, effective mass megp, and peak amplitude x) = 
J 2kgT / Met Q?, at temperature 7. Integration yields 
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+092, 2 


XO 
ay e Tt/2 


‘m 


(17) 


sinQmt, 


where we have used the high-Q approximation (Qn >> IT) to 
neglect the cosine term. 
Substituting Equation into Equation (15), we have 


ay (t) © aye (3-4)! (1 + o2or [ein = ent) , (18) 


where we have written 


_ Gxo 
c= om (19) 


and used the Jacobi-Anger expression to write 


. : pkey Lj 
eiAsing _ 3 Ji (A)elk? ~1+iAsing (20) 


k=—00 


for small A = €e"/? < 1. Here, J;(A) are the kth Bessel 
functions of the first kind. 

We justify this approximation by first noting that e~™ 2] 
for all finite positive times and then rewriting € as 


kg 
280\/ fim _ 2g0vi 


= Qn Qn 


(21) 


using the vacuum optomechanical coupling rate go = Gxgp¢ 
and the amplitude of the mechanical zero-point fluctuations 
Xopt = Vh/2mep~Qm. Here, i = kgT/hQm is the average 
phonon occupation of the resonator for kgT >> hQm. In this 
form, go < Qm is absolutely necessary to have € < 1, and 
we must additionally consider the phonon occupation 7. For 
optical frequency optomechanical devices, this condition is 
commonly satisfied, especially at low temperatures where 7 
is small. For the bottle resonators used in our experiments, 
go /2m ~ 350 Hz, so € <1 and thus Ee!"/? < 1 as long as 
i < 10!° (corresponding to a temperature of 10’ K). 

We now return to Equation to look for the particular 
solution to Equation (ip. We note that 


G(t) = an (t) +an(t) F(t) +an() FO 


=~ (F —iA—iGx(0)) (an(t) +ay()F() + V Resin): 
(22) 


Given that the homogeneous solution obeys Equation (14), it 
follows that 


_ Ves) 


P(t (23) 
me) a(t) 
For the mechanical mode thermometry, we phase modulate 
the input laser light by driving an electro-optic modulator with 
a sinusoidal signal of the form Yor "/ 2 sin Qmoat, where Vo is 


the drive voltage amplitude, Qyoq is the drive frequency and 


y < Tis the linewidth of the driving source. This produces an 
input field of the form 


iBe~Vt/2 si 
Sin (t) = Sine PS sinQ oat 


oo ; 
=a) Fp(Be Peet 


k=—0o 


(24) 


where we define the phase modulation depth as B = 2Vo/Vz, 
given the device-dependent half-wave voltage, V;. If the 
phase modulation is weak (6 < 1), it is sufficient to consider 
only the first order sidebands at OQmoq. In this case, we can 
again use the approximation in Equation to write 


sin(t) = Sin (1 + Bern -  _iSneat = eos ) . (25) 


We substitute this result, along with the homogeneous so- 
lution of Equation (18), into Equation (23), yielding 


f(t) _ V KexSin ,(§—iA)t (1 + Bem ere = cont) 


ao 


. (1 _ ge [ein ~<a) | 


Keeping only terms to first order in the small parameters € and 
B, integrating Equation and multiplying by ay(t) (Equa- 
tion (18)), we obtain the particular solution 


(26) 


B yt e!Qmoat 
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: e2 ; 
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eQmt eiQmt ef Qmt 
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(27) 


where we have used the fact that k >> IT, y to simplify the final 
expression. We note that a;,(t) decays much more rapidly than 
dp(t), so we neglect ap(t) and take a(t) ~ ap(t). 

We now calculate the field output by the cavity by substi- 


tuting Equations and into Equation (12), 


, ‘ K 
Sout(t) =Sin (1 + Ben [einen a aa _ ais 


7 
Kex B Yt | ei Qmodt eQmoat | 
_ nade 


2 _i(A—Qmoa) % —i(A+Qmod) 
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e 


(28) 


B. Detection 
In our experiments, we use a direct detection scheme where 


a photodetector sensitive to the intensity of the light outputs a 
voltage signal V(t) proportional to |sout(t)|*, 


V(t) = Voc + Vin(t) + Vinoa (t). (29) 


We have separated the output signal into a part 


Voc = H(0)rin (1- a (30) 


2 2 
ee 


which is constant in time, and two high-frequency parts which 
oscillate at the mechanical resonance frequency, 


= 5. § -rr/2 Ci (Qm,A) cos Qmt 
Vin{t) $2 Om)fia ker € "| D0,A)D(On,A)D(—On, A) 
a C2(Qm, A) sinQunt 
D(0,A)D(Qm,A)D(—Qm, A) )’ 


(31) 
and at the modulation frequency, 
= . B —yt/2 
Vinod (t) = AH (Qmod)Min Kex g° x 
( Ci (Qmoa, A) cos Quoat (32) 
D(0,A)D(Qmoa, A)D(—Qmoa; A) 


' Cr (Quod , A) sin Qunoat ) 
D(0, A)D(Qmod)D(—Qmod) 


respectively. H(@) is a function which describes the fre- 
quency response of the photodetector and other detection elec- 
tronics. We assume that it is locally flat, such that for a suit- 
ably chosen Qmyog, we can take H (Qmoa) = H(Qm) =H. Note 
that Vpc is simply the transmission profile for the optical reso- 
nance, as given by Equations (1) and (2) in the main text. For 
convenience, we have also defined 


C,(Q,A) = AQ(—K[K? +4A7] + Kx [K? +4A? —407]), (33) 


C)(Q,A) = AQ? (K[4K., — 3x] +4[A? — 7), (34) 


and 


K2 
D(Q,A) = > + (A—Q)?. (35) 


To obtain the power spectral density of the voltage signal, 
we must first Fourier transform V(t). In our experiment, this is 
done as part of the software post-processing, but it can also be 
performed in hardware, using a network or spectrum analyzer, 


for example. Vpc is filtered out by a high-pass filter on our 
photodetector so we address only the high-frequency part of 
V(t). 

We define the Fourier transform as 


co 


V(@) = / eV (t)dt. (36) 


—oo 


Since the Fourier transform is linear, we apply it to V(t) and 
Vinoa(t) separately. We begin by substituting the form of & 
given by Equation into Equation and noting that 


LT. 
X(t) = — = cos Ont — Ome sin Qt (37) 
—It/2 


= —Omx0e sinQ pf. 


Equation then becomes 


Mi = (S (Bardi) ClBn. AC /9n) | 


roe D(0,A)D(Qm,A)D(—Qq, A) 
(38) 
with its Fourier transform given by 


Vin(@) = Any 


Gkex (C1 (Qm,A) — i@C>(Qm,A)/Qm 
22m 


Near the mechanical resonance frequency, @ © Q, and 


Vin(@) = Anin 


GKex ( Ci (Qm, A) — iC2(Qm,A) )x(o) 
20m \D(0,A)D(Qm,A)D(—Qm,A) 
(40) 
Similarly, Vinoa(t) can be rewritten in terms of the original 
phase modulation signal @(t) = Be~”/? sin Qmoat, 


Hn: Kex ( C (Qmoa, A) (t)/Qmoa 
" 2 D(0,A)D(Qmoas A)D(—Qmoas A) 
a C2 (Qmod; A)o (t) ) 
D(0,A)D(Qmnoa; A)D(—Qnoa; A) 


Vmod (t) = 


(41) 


Analogous to Equation (40), the Fourier transform of Equa- 
tion (41) near the modulation frequency is 


K, iC} (Qioa , A) 


Vinod (@) = Anin 5) me Sa 
C2(Qmoa, A) 
fo) 
ee See (), 


(42) 


where ®(@) is the Fourier transform of @ (ft). 
Combining the results of Equations and (42), we cal- 
culate the voltage power spectral density, 


D(0,A)D(Qm,A)D(—Qm.A) ) x(a), 


= 
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FIG. 8: Calculated ratio of @(Qm,A) to @(Qmoa, A) for experimental 
parameters (K/27 =50 MHz, kex /27 = 20 MHz, Qn /2a =55 MHz 
and Qmoa/27 = 54.5 MHz). A variation of less than 3% in this ratio 
across the optical resonance is shown. 


. il 
Syv(@) = lim —|Vn(@) +Vmoa(@)|”>; (43) 


which can be reduced to 


wt bc i 
Svv(@) = lim =|Yn(@)|? + lim —|Vmoa(@)|? (44) 


if we choose Qmoq such that |Qmog — Qn| > T,y. This pre- 
vents overlap between the peaks in the spectrum and allows 
us to neglect the cross terms. The final result is then 


G2 
Syv(@) = ne, (= 2 (Qm, A)Sxx(@) + &(Qmoa, A)S oo ()) 
7 (45) 
where S49 (@) is the PSD of the phase fluctuations induced by 
the phase calibration signal. We have defined the transduction 
coefficient as 


H?Ke —-C7(Q,A) +C3(Q,A) 
4 D?(0,A)D2(Q,A)D2(—, A)’ 


a:(Q,A) = (46) 


We see that for a choice of the modulation frequency 
Quod & Qm, O(Qmod;A) ¥ O(Qm,A) and the phase modula- 
tion signal is transduced nearly identically to the mechanical 
motion by the optical cavity and photodetector. In actuality, 
a is also dependent on the detuning of the laser, but it can 
be shown that for |Qm— Qmoa|/Qqnm ~ 1%, the variation of 
the ratio of &(Qm,A) /O(Qmoa, A) with detuning is negligibly 
small. For the parameters in our experiment, Figure [8] illus- 
trates that this variation across laser detuning is less than 3% 
in our experiment. 


Although we have focussed on direct detection techniques, 
this method of calculating the detected PSD applies equally 
well to other detection schemes, including optical homodyne 
or heterodyne systems. It can be shown that analogous results, 
albeit with a different functional dependence on the laser de- 
tuning, can be obtained for such systems. 

We finally note that we have made no assumptions about 
the optomechanical system, beyond requiring that k >> T, as 
is true for any system in the standard hierarchy of optome- 
chanics. The only constraints placed on the modulation sig- 
nal are that it is small (8 < 1) and that its frequency is cho- 
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